DESY 95-162 



Measurement of the Total Photon Proton Cross Section 

and its Decomposition 
at 200 GeV Centre of Mass Energy 

HI Collaboration 



Abstract: 

We present a new measurement of the total photoproduction cross section per- 
formed with the HI detector at HERA. For an average centre of mass energy of 
200 GeV a value of a^ t = 165 ± 2 ± has been obtained. A detailed analysis of 
the data in adequate kinematic regions enabled a decomposition of the total cross sec- 
tion in its elastic, single diffractive dissociation and remaining non-diffractive parts, 
based on safe assumptions on the double diffractive dissociation contribution. 
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1 Introduction 



The total cross section is an important quantity related to the fundamental properties of particle 
interactions. Although measurements are available for hadron hadron and real photon hadron 
collisions at low energy , so far only the data from pp-colliders provide precise information on 
the rise of the total cross section at high energy. The ep collider HERA, with e and p energies of 
27.6 and 820 GeV, provides a new source of information on high energy photon proton collisions. 

The interaction of electrons and protons at the HERA collider is dominated by photoproduc- 
tion processes, in which the electron scatters through small angles emitting a quasi-real photon, 
which then interacts with the proton. Recently the total photoproduction cross section has been 
measured at HERA at a jp CMS energy of 195 GeV by HI § a™ = 156 ± 18 fib Q and at a jp 
CMS energy of 180 GeV by ZEUS || aj^ t = 143 ± 17 fib. These measurements confirmed the 
expected rise of the total "yp cross section with energy. However, large systematic uncertainties 
do not yet allow discrimination between different models predicting a moderate rise of the "yp 
total cross section [4-8]. In the previous HI analysis the systematic error is dominated by 
the assumptions on the partial 7p cross sections, which is important as they have different ac- 
ceptance. ZEUS H determined the fraction of -yp diffractive events directly from the data and 
thus reduced the model dependence of the result. However, the precision of this measurement 
is limited by the large systematic error in the efficiency of tagging the scattered electrons. 

In this paper a new determination of cr/o* at the average centre of mass energy W 1V = 200 GeV 
is presented. Dedicated trigger conditions, specially designed for high energy photoproduction, 
allow the measurement of the diffractive components of the "yp cross sections and, therefore, 
substantially reduce the model dependence in the acceptance calculations. This, together with 
a better understanding of the electron tagging efficiency improves the accuracy of the cr/ot mea- 
surement compared with the earlier results from HERA. 

The data used in the present analysis were taken during a period in which HERA was 
operated with a positron beam. Nevertheless, "electron" is used as generic term for the HERA 
beam lepton throughout this paper. 



2 Photoproduction 

In ep collisions the total photoproduction cross section, aj^ t can be related to the total differential 
ep cross section by the Weizsacker- Williams formula |J for the photon flux F(y,Q 2 ) 

r l2 ep( \ 

aZ(ys)-(l + S RC )-F(y,Q 2 ) = 



dydQ* 



2irQ?\ y y Q J 



where Q 2 is the negative square of the photon 4-momentum, or the virtuality of the photon, 
and s is the squared centre of mass energy of the ep interaction. For small scattering angles, y 
is defined as 1 — E' e /E e , where E e and E' e are the energies of the initial and scattered electron 
respectively. The minimum photon virtuality is Q min = (m e y) 2 /(l — y). The factor (1 + 5rc) 
takes into account QED radiative corrections to the ep Born cross section. In the formula (|l|) a 
dependence of u^ot on Q 2 and a contribution of longitudinally polarized photons are neglected. 
These are good approximations in the present kinematic conditions |T~ 



Note that this number underestimates aj^ t by 6-7% due to the approximation used in the theoretical expres- 
sion for the photon flux. 



Hadronic final states produced in real photon proton collisions resemble those observed 
in hadron hadron collisions. This similarity led to the phenomenological approach to describe 
photoproduction by the so called Vector Meson Dominance (VMD) model [JO], where the photon 
first converts into a vector meson (predominantly the p°) which then interacts with the proton. 
As in hadron hadron collisions, the total 7p cross section has a substantial contribution from 
diffractive 7p reactions, which have a final state topology radically different from the bulk 
of non-diffractive events. Since diffractive reactions involve no exchange of quantum numbers 
between the incident particles, the final state is characterized by the appearance of large rapidity 
intervals, or gaps, with no hadrons. This feature of diffractive events is exploited below to 
determine their contributions to the jp cross section. We distinguish the following diffractive 
processes in photoproduction: 



• Elastic vector meson production (EL) j + p — > V + p, where V stands for one of the vector 
mesons p°,ui,(fi. The true electro-magnetic elastic reaction 7 + p — > 7 + p has a very low 
cross section and is neglected [11]. 

• Single photon diffractive dissociation (GD) 7 + p — > X + p, where the photon dissociates 
into the heavy hadronic state X and the proton stays intact. 

• Single proton diffractive dissociation (PD) 7 + p —* V + Y, where the proton dissociates 
into a hadronic state Y and a vector meson is produced in the photon direction. 

• Double diffractive dissociation (DD) 7 + p — > X + Y, where both the photon and the 
proton dissociate. 



For the last three diffractive reactions the cross section is considered for the full momentum 
transfer range and for masses (M) of the dissociating system obeying M 2 < 0.1W 2 p , where 
Wryp is the centre of mass energy. The value 0.1 is chosen to enable direct comparisons with 
measurements of the diffractive cross section from hadron colliders and fixed target photopro- 
duction experiments as well as with available theoretical calculations. 

All processes 7 + p — » X not belonging to the contributions defined above are called non- 
diffractive (ND). These processes dominantly involve exchange of quantum numbers between 
the photon and the proton. 



3 Monte Carlo Models for Photoproduction 



Two Monte Carlo (MC) models, based on the event generators PYTHIA |2j and PHOJET [jj, 
are used for the acceptance calculation. Both models include all the diffractive and non- 
diffractive contributions to the jp cross section discussed above. 

The model for non-diffractive events in the PYTHIA program is similar to the multiple- 
interaction model developed for hadron hadron collisions O]. An eikonal approach is used, in 
which the rate of jets above a transverse momentum p± m i n (the default value p_i_ m m = 1-45 GeV/c 
is taken) is combined with a parameterization of the non-diffractive total cross section to cal- 
culate a probability distribution in the number of semi-hard interactions. For events below this 
cut-off two longitudinal strings are stretched between the proton and the "VMD-photon" , to 
give a representation of an event structure caused by soft gluon exchange. 

Events with elastic, diffractive single and double dissociation scattering in the PYTHIA MC 
model have the same general structure: the t-dependence is given by a function exp(Bt), where 
t is the square of the four-momentum transfer in the diffractive reaction and B is the nuclear 
slope parameter. Within the energy range used in this analysis the nuclear slope for elastic 



vector meson production is B ~ 11 (GeV/c) -2 . A p° formed by 7 — > p° in elastic or diffractive 
scattering is transversely polarized and therefore its decay angular distribution in p° — > 7r + 7r~ 
is taken to be proportional to sin 2 9, where the reference axis is given by the p° direction of 
motion. The relative rates of p°,uj,4> production are assumed to be about 13 : 1.5 : 1 [11]. 



In single diffractive dissociation, the -B-slope is assumed to be half that of elastic scattering. 
The events are generated according to a dM 2 /M 2 distribution for the dissociation system of 
mass M, The mass spectrum of the system is assumed to start at 0.2 GeV/c 2 above the 
mass Mi n of the incoming particle (using the p° mass for the incoming 7). A light dissociated 
system, with a mass less than 1 GeV above the mass of the incoming particle, is taken to 
decay isotropically into a two-body state. Single-resonance states, such as N* or w(1600), are 
not generated explicitly, but are described in this average manner. A more massive system is 
treated as a string stretched along the 7^ interaction axis. The secondary hadrons from the 
string decay are distributed in a longitudinal phase space with limited transverse momentum. 

In the event generator PHO JET, the multi-particle non-diffractive final states are constructed 
from a parameterization of the photon proton scattering amplitude in an eikonal approximation 
using the two-component Dual Parton Model [19]. The coding of the model is similar to that 



of the MC generator DTUJET [16] simulating particle production in pp and pp collisions up to 
very high energies. 

In the generator PHOJET elastic vector meson production is similar to that in the PYTHIA 
model. For diffractive dissociation the PHOJET model assumes a mass dependent nuclear 



slope B [17|. This slope parameterization gives a steady transition from elastic scattering 
to single and double diffractive dissociation. The mass spectrum is generated according to 
a dM 2 /(M 2 — M 2 n ) distribution starting from two pion masses above the mass Mj n (using the 
p° mass for the incoming 7). The low-mass resonance structure is taken into account in an 
approximate way to provide a phenomenological description of the dissociated mass spectrum 
observed experimentally |l8|, O]. To take the transverse polarization of the incoming photon 
into account, the decay of the elastically produced vector meson resonances into two or three 
particles is performed in the s-channel helicity frame according to the angular distributions 
given in [pi]] . In addition to resonances, a continuous multi-particle final state in diffraction is 
generated by simulating a pomeron-proton or pomeron-photon scattering exactly as in the Dual 
Parton Model used for photon hadron scattering. The pomeron is treated like a virtual meson. 
The soft and hard scatterings in diffraction are generated according to cross sections given by 
Regge-parameterizations and the QCD Parton Model, respectively. 

The generated events are fed into the HI detector simulation program and are subject to 
the same reconstruction and analysis chain as the real data. 



The QED radiative corrections are calculated using the HERACLES MC program |2(J . This 
takes into account single photon emission from the lepton line as well as the self energy correction 
to the Born photoproduction cross section. 



4 HI Detector 



A detailed description of the HI apparatus can be found elsewhere [21]. A schematic layout of 
the central HI detector components is shown in Fig. |]. In the following we briefly describe the 
components of the detector relevant for this analysis. 



Measurements of charged particle tracks and the interaction vertex are provided by cen- 
tral and forward tracking systems, both consisting of drift and multi-wire proportional cham- 
bers (MWPC). The central and forward track chambers cover the complete azimuthal range and 
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Figure 1: A layout of the central part of the HI detector. 



—2.0 < i] < 3.0 in pseudo-rapidity rj = — ln(tan |). Here 9 is the polar angle with respect to the 
proton beam direction (positive z axis). The central jet chamber (CJC) is interleaved with inner 
and outer double layers of MWPC, which were used in the trigger to select events with charged 
tracks pointing to the interaction region. This MWPC system covers the range —1.5 < rj < 1.5. 
A backward proportional chamber (BPC), with an acceptance of —3.0 < r] < —1.5 allows efficient 
detection of charged particles produced at large 9. 

The tracking region is surrounded by a fine grained liquid argon (LAr) calorimeter 22] 
consisting of an electro-magnetic and a hadronic section. The total depth of the LAr calorimeter 
varies between 4.5 and 8 hadronic interaction lengths. Under test beam conditions it has an 
energy resolution cr/E ra 12%/y/E/GeV © 0.01 for electrons and » 50%/y/E/GeV © 0.02 for 
pions. The LAr calorimeter covers the complete azimuthal range and —1.5 < rj < 3.3. The 
calorimeter is surrounded by a super-conducting solenoid providing a uniform magnetic field of 
1.15 T parallel to the beam axis in the tracking region. 

The time of flight system (ToF) is located at z ~ —2m behind the Backward Electro- 
Magnetic Calorimeter (BEMC), which is about one hadronic interaction length deep. ToF is a 
hodoscope consisting of two planes of plastic scintillators mounted perpendicular to the beam 
direction. The angular coverage of the ToF counters corresponds to —3.5 < r\ < —2. Having 
a time resolution better than 2 ns the ToF system enables efficient separation of ep interaction 
events from the upstream background. In the present analysis the ToF system is included in the 
trigger for photoproduction events. The efficiency of the ToF counters has been measured using 
muons in the proton beam halo and corresponds to (98 ± 1)% for minimum ionizing particles. 

The luminosity system measuring the reaction ep — > e^p consists of two TICl/TlBr crystal 
calorimeters. The small angle electron detector (electron tagger) is located at z = —33 m and is 
also used to trigger on photoproduction events. Its 7 x 7 crystal matrix (an individual crystal 
measures 2.2 x 2.2cm) accepts electrons with an energy between 0.2E e and 0.8-Ee and scattering 
angles 9' < 5mrad (9' = tt — 9), corresponding to Q 2 < Qmax = 0.01 GeV 2 . The photon 
detector is located at z = —103 m and consists of a 5 x 5 crystal matrix. Both calorimeters 
are 22 radiation length deep and their energy resolution in the present data taking period was 
measured to be a{E)/E = 0.15/yE/GeV © 0.01. 



5 Luminosity Measurement and Tagging Efficiency 



The basic requirement used to tag quasi-real photoproduction processes in HI is the detection of 
the scattered electron in the electron tagger. This guarantees very low Q 2 < 10~ 2 GeV 2 . Three 
main ingredients contribute to the overall precision of the measurement of the cross section crj^\ 
the luminosity measurement error, the knowledge of the electron tagger acceptance and the 
efficiency of the main apparatus for triggering and reconstruction of the hadronic final states 
produced in photon proton collisions. The first two are discussed in this section. 

The luminosity measurement utilizes the Bethe-Heitler (BH) process ep — > cyp [^3]. Several 
methods for the measurement of the BH process can be exploited and they are described in detail 
in |24j. In this analysis the method based on the measurement of the photon energy spectrum 
with E~ > E m i n ~ 8 GeV is used. The value of E m i n is chosen such that it is well above the 
photon detector trigger threshold. The error in the luminosity measurement is then dominated 
by the precision of the energy calibration and by the correction for the complex structure of 
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Figure 2: Performance of the HI luminosity system and the electron tagger acceptance, (a) 
Normalized photon energy distribution for Bethe-Heitler events in the data (symbols) and Monte 
Carlo with (full line) and without (dashed line) the event pile- up effect, (b) Correlation between 
the energy and the lateral coordinate of the impact point of scattered electrons in the electron 
tagger in the data (symbols) and Monte Carlo (full line). Vertical lines indicate the fiducial 
cut \xet\ < 6.5 cm used in the analysis, (c) Energy spectrum in the electron tagger for Bethe- 
Heitler events; data (symbols) are compared to MC with measured e-beam tilt X = —0.13 mrad 
and different offset values of the electron trajectory in the HI interaction point: x Q g = -0.5mm, 
+0.5mm, +1.5mm (dashed, full and dashed-dotted lines respectively), (d) y-distribution in the 
electron tagger for jp events in the range 0.3 < y < 0.7 used in this analysis. The points 
represent the HI data, histograms are Monte Carlo predictions for the models PYTHIA (solid) 
and PHOJET (dashed). 



the proton bunches - the so called "satellite bunch" effect p4|] . After the final absolute energy 
calibration of the luminosity detectors the observed photon energy spectrum is shown in Fig. |||a. 
It is well described by the BH process simulation taking into account energy resolution and pile- 
up effects (ie. several overlapping ep — ► cyp events in the same bunch crossing). The precision 
of the integrated luminosity measurement in different 1994 data samples varies between 1.5% 
and 5.6%. 

The electron tagger performs a double function. It both tags 7p events and provides a 
measurement of the scaling variable y. Since the precision of the energy and coordinate recon- 
struction is not good enough in the areas close to the detector boundaries, xet = ±7.7 cm, a 
fiducial cut \xet\ < 6.5 cm is used in the analysis (Fig. ^b). A corresponding cut on \vet\ is 
redundant because of the confinement to the HERA bending plane. 

The acceptance A(y,Q 2 ) of the electron tagger for scattered electrons depends strongly 
on the HERA electron beam optics, being most sensitive to the horizontal tilt 6 X and the 
horizontal offset Ax ff with respect to the reference trajectory at the HI interaction point. The 
electron beam tilt (typically 0.1 mrad) can be measured with a precision of ±0.02 mrad by 
monitoring the position of the photon spot at the photon detector. The offset is not measured 
directly. The acceptance A BH (y) , integrated over Q 2 , can be determined using ep — > ejp events. 
However, A BH (y) differs from the acceptance for photoproduction A(y) due to the different Q 2 
dependences. 

The following procedure has been used to determine A{y) for any data sample with constant 
beam conditions. First, the acceptance A BH (y) was measured from ep — ► ejp events. Then the 
Monte Carlo program simulating the HI luminosity system together with the HERA beam optics 
was tuned to the data by varying Ax g, which is the only free parameter in the procedure. Fig. ||c 
illustrates the sensitivity of the energy distribution in the electron tagger (and thus A (y)) 
to the horizontal offset of the reference trajectory. A precision of Ax g = ±0.2 mm has been 
achieved by this procedure. Finally, the acceptance A(y) was calculated using the measured tilt 
6 X and the tuned value of Ax fr- The errors have been estimated from calculations of A BH \y) 
using extreme values of the parameters. A limited range of 0.3 < y < 0.7 was used in the 
analysis, to avoid tails where the acceptance value is less than 20%. Within this interval, errors 
between 3% and 5% were obtained in the value of J A{y)dy for different data samples. We 
therefore conclude that 5% can be used as a conservative estimate of the precision to which the 
electron tagger acceptance is known in this analysis. Fig. ||d shows the comparison of the y 
distributions in the data with Monte Carlo, using the two different models for photoproduction 
as described in section 3. 

6 Trigger Conditions and Event Selection 

The data used for the measurement of a]^ t were collected in a short dedicated period during the 
1994 data taking. The HERA machine was operated with 153 colliding bunches of 27.6 GeV 
positrons and 820 GeV protons. In addition, 32 "pilot" bunches, 17 proton and 15 positron, 
had no counterpart and produced no ep collisions enabling an estimate of the beam induced 
background. Two data samples were collected with different mean z-positions of the ep interac- 
tion vertex: the nominal position at z = 4 cm and a position shifted in the proton direction at 
z = 71 cm (see Fig. [J). 

The data collected with the nominal vertex position correspond to the integrated luminosity 
of 23.8 ± 0.4 nb _1 , while the shifted vertex position data correspond to 23.8 ± 1.3 nb -1 . The 
advantage of the shifted vertex data is the higher acceptance for diffractive reactions in the 
region where diffractive processes can be separated safely from the majority of non-diffractive 
events. 



The data were taken with two independent trigger conditions. The first trigger condi- 
tion, termed "ToF-trigger" , is formed by the coincidence of a signal from the electron tagger 
(E' e > 4 GeV) with a signal from the ToF system coming within the time interval expected 
for ep interactions. The ToF-trigger is fired by hadrons originating from photon fragmentation. 
This trigger is efficient for all classes of photoproduction events including elastic vector meson 
production, although it is affected by the BEMC material in front of the ToF system. The ToF 
trigger was enabled in both data samples. 

The second trigger condition, termed "Ray-trigger", requires a coincidence of the electron 
detector signal with at least one track pointing to the vertex region. The track condition is 
derived from the cylindrical MWPC and requires &pt ~ 200 MeV/c. This trigger has been used 
in previous HI analyses [§]. In the present analysis the Ray-trigger is used for cross checks and 
was activated during the run with nominal vertex position only. 

Both Ray- and ToF-triggers require in addition the energy in the photon detector to be less 
than 2 GeV. This condition substantially reduces the size of QED radiative corrections and also 
suppresses accidental coincidences of Bethe-Heitler events with beam induced background. 

The triggered events are subjected to several offline cuts. The fractional energy of the 
photon, as measured by the electron tagger, is required to be in the interval 0.3 < y < 0.7. 
The event vertex, reconstructed from tracks in the CJC, must be within ±30 cm of the mean 
z-position of the interaction point. The vertex z-position distribution has a Gaussian shape 
with a sigma of 10 cm reflecting the length of the proton bunches. For ToF-triggered events at 
least one reconstructed track in the CJC or in the BPC is required in addition. Similarly the 
Ray-triggered events are required to have at least one CJC track in the region covered by the 
Ray-trigger. 

Several sources of background contribute to the data samples. The main background source 
in the ToF-triggered event samples is electron interactions with residual gas ("beam-gas") or 
with material inside the beam-pipe ( "beam- wall" ) . This contribution is estimated using the data 
from the non-colliding (pilot) electron bunches and amounts to 4% and 8% respectively, in the 
event samples with shifted and non-shifted vertex. For Ray-triggered events this background is 
negligible. Another major source of background originates from accidental coincidences of the 
electron tagger signal with events resulting from proton beam-gas collisions within the nom- 
inal ep interaction region. This contribution is estimated using special monitoring triggers, 
with looser triggering conditions, and is about 3% in the Ray-triggered data sample. For the 
ToF-triggered events this background is negligible. Still another type of background stems from 
the accidental coincidence of electron beam induced background with a Bethe-Heitler process in- 
duced signal in the electron tagger. This coincidence appears as background when the associated 
photon escapes detection in the photon detector, due to the small inefficiency of this detector. 
From the measured rate of the Bethe-Heitler process this background contribution is estimated 
to be about 1% in all data samples. Finally there is a small background contribution in all data 
samples from the QED 2-photon lepton pair production processes, with one photon emitted 
from the incident electron and the other photon emitted from the proton. This background was 
calculated using the LPAIR MC event generator |25| and amounts to less than 0.2% under the 
present trigger and selection conditions. The selected event samples vary between about 19,000 
and 22,000 events. The background is subtracted on a statistical basis in the analysis. 

Fig. [3| shows a comparison between the nominal vertex ToF-triggered data sample and the 
two MC simulations. In Fig. ||a the z-position of the reconstructed event vertex is shown and in 
Fig. |3]b-d the multiplicity, transverse momentum and pseudo-rapidity of charged tracks are given. 
Events with no CJC tracks (i.e. events which fulfil the BPC requirement) do not contribute in 
the distributions of Fig. ||a, c and d. The agreement between data and simulations demonstrates 
that the MC event generators reproduce well the main features seen in the data. Especially the 
z-vertex distribution in Fig. |||a shows that the background contamination in the data is small. 
Similar comparisons (not shown) for the Ray-triggered event sample lead to the same conclusion. 




Figure 3: The measured z-position of the reconstructed event vertex (a), multiplicity (b), 
Pt (c) and rj (d) distributions of charged tracks in the data (points) for the nominal vertex 
sample compared with those in MC simulation using PHOJET (full histogram) and PYTHIA 
(dotted histogram). 

7 Cross Section Calculation Method 

For each data sample the observed number of events N is related to the differential ep cross 
section (|l|) by the expression 

dN r r ^ a, ^ d 2 a e P , , 
C-e{y)- A{y,Q 2 ) (2) 



dydQ 2 Ka ' K ° 7 ^ ' dydQ 2 

where C is the integrated luminosity, e(y) is the efficiency of the trigger and selection criteria 
for the main HI detector, and A(y, Q 2 ) is the acceptance of electron tagging as described above. 
The A(y, Q 2 ) is factorized out since it depends only on the parameters of the scattered electron 
and not on the details of the photoproduction process. We assume e{y) to be independent 
of Q 2 since the transverse momentum of the scattered electron is always small. Integrating 
equation (||) over y and Q 2 in the range from y m in to y ma x and from Q 2 min to Q 2 max gives 

N=(l + 6 RC )-£-F-A-E-aZ, (3) 

using the photon flux integral F and averaged values ^ of A, e and the cross section a^ t . An 
identical expression also holds for any partial cross section dj with only e depending on the 

2 A strict definition of the quantities averaged over y and Q 2 is 
<t = I crZ(ys)e(y)A(y)F(y)dy/ J S (y)A(y)F(y)dy, F = J F(y)dy, A = J A{y)F(y)dy/F, 
e = J e(y)A(y)F(y)dy/AF, with F{y) = J F(y,Q 2 )dQ 2 and A(y) = / A(y, Q 2 )F(y, Q 2 )dQ 2 /F(y) 



sub-process i. Therefore, one has for their sum 



N = (l + 8 RC )-£-F-A-J2^i- (4) 

Eq.(|4|) is also valid for the number of events in any kinematic domain of the photoproduction 
process, with e, determined accordingly. This enables us to find the partial cross sections by 
considering appropriate kinematic regions, enriched by different sub-processes, and solving a 
(generally over-constrained) system of equations for the Oi , with efficiencies e$ calculated by MC 
simulations. 

For our basic kinematic limits of y m in = 0.3 , y ma x = 0.7 (corresponging to the range 
165 < W 1P < 252 GeV) and Q 2 max = 0.01 GeV 2 we find F = 0.0136 and the acceptance A 
varying from 0.546 to 0.570 depending on beam conditions. This variation of A is properly 
taken into account in the analysis. The MC efficiencies £{ for the selected ToF-triggered samples 
are presented in Table [j] . The average value of the radiative correction 6rc is estimated to be 
(— 1 ± 1)% and (+1 ± 1)% for the ToF- and Ray-trigger samples, respectively. 

Table 1: Efficiencies £{ (%) for the different subprocesses of jp scattering as calculated using 
the Monte Carlo simulation based on the PHOJET(PYTHIA) models for various data samples. 
The rj max < and r] m i n > 1 samples are used to find the diffractive contributions. 



sample 


sub-sample 


subprocess 


GD 


PD 


DD 


EL 


ND 


shifted vertex 
ToF trigger 


all events 

Vniax *v 


74(77) 
31(28) 
0.2(0.3) 


68(71) 
28(18) 
21(26) 


71(73) 
15(17) 
3.4(3.1) 


57(54) 
52(50) 
0.4(0.4) 


61(65) 
0.1(0.1) 
0.1(0.1) 


nominal vertex 
ToF trigger 


all events 


66(73) 


42(46) 


68(72) 


29(26) 


65(70) 


nominal vertex 
Ray trigger 


all events 


57(65) 


8(14) 


62(52) 


0(0) 


95(94) 



8 Cross Section Measurement 

To measure the total 7p cross section we use data with the nominal vertex position, where the 
uncertainty in the luminosity calculation is significantly smaller than in the data taken with the 
shifted vertex position. However, as can be seen from Table [l] the efficiencies for the diffractive 
channels are higher in the shifted vertex data sample. Therefore, for the determination of the 
diffractive contributions to the 7p cross section we use the data taken with shifted ^-vertex 
position. The analysis and the cross checks are described in the following subsections 8.1 — 8.3. 

8.1 Diffractive Contributions 

To measure the diffractive 7p cross sections we choose the variables in which a separation of 
diffractive and non-diffractive contributions is least model dependent. These variables, rjmax and 
r/min, are related to the central rapidity gap in the hadronic final state. 

For each event r] max is defined as the maximum pseudo-rapidity of all reconstructed charged 
tracks and all clusters in the LAr calorimeter with energy larger than 400 MeV. In diffractive 



events rj max indicates the maximum pseudo-rapidity of secondary hadrons from photon dissocia- 
tion. It was shown in a recent HI analysis of diffractive photoproduction (2(| that the spectrum 
of 

Vmax for non-diffractive events falls nearly exponentially with decreasing rj max , whilst the rate 
of photon diffractive dissociation depends only weakly on this variable. Events with elastic vec- 
tor meson production have the largest possible width of the rapidity gap and are concentrated 
at lower values of the rj max spectrum. 
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Figure 4: The rj max (a) and rj m i n (b) distributions for the sample with shifted vertex. The 
data are shown with solid circles. The histogram represents the result of a combined fit based 
on the PHOJET model summing all five contributions (EL, PD, GD, DD and ND) and using 
udd = 20 /ib. The full curve in (a) represents the sum of elastic and single proton dissociation. 
The sum of single photon and double dissociation contributions in (a) is represented by the 
dashed curve. The sum of single proton and double dissociation in (b) is shown by the full 
curve. 

The 

Vmax distribution for the data sample with shifted interaction vertex is shown in Fig. |4|a, 
where we additionally require the calorimetric energy with rj > 1 to be less than 1 GeV. This 
requirement reduces the contribution from non-diffractive events. MC calculations show that all 
four diffractive reactions contribute to the rj max spectrum in the range —3.5 < r] max < 0, whilst 
the contribution from non-diffractive processes is negligible. Below r] max = —2 the spectrum 
is dominated by elastic and single proton dissociation channels, where rj max is determined by 
the maximum pseudo-rapidity of the vector meson decay products. These two contributions 
are practically indistinguishable by shape, but have different efficiencies. For proton diffrac- 
tive dissociation only the events with a low mass proton system (M < 10 GeV/c 2 ), where 
secondary particles escape detection very close to the proton beam direction, contribute to the 
Vmax spectrum. The region —2 < ri max < is dominated by single photon and double diffractive 
dissociation contributions, which have again a similar shape, but different acceptances. 

The 

ijmin variable is defined as the minimum pseudo-rapidity in the interval —2 < rj < 3.5, 
of all charged tracks and of all calorimeter clusters with an energy larger than 400 MeV. An 
additional condition is that the event has a reconstructed charged track with rj < —2.4. The 



latter condition is necessary for tagging a low mass hadronic system on the photon side. In 
diffractive events ry m j n indicates the minimum pseudo-rapidity of secondary hadrons from proton 
dissociation. The measured r] m i n distribution for the data sample with shifted interaction vertex 
is shown in Fig. |]b and is dominated by the non-diffractive contribution, which falls nearly 
exponentially with increasing r) m %n- However, there is a flat part of the spectrum with T) m % n > 1 
caused by proton diffraction dissociation processes, where the mass of the hadronic system 
produced on the proton side is larger than about 5 GeV/c 2 . Both single and double dissociation 
contributions have a similar shape. The acceptances for various partial processes to contribute 
to the r] max and r\ m i n distributions are presented in Table |], calculated using the PHOJET and 
PYTHIA MC models. 



hum 



To obtain the diffractive 7p cross sections wg make a, combined fit of the t]rnax 

and r} Tl 

distributions in the intervals —3.5 < r] max < and 1 < r/ m i n < 3.5 using formula @). In this fit 
procedure the shapes of the spectra for each partial contribution are fixed from the MC calcu- 
lations, while the cross sections U{ for single dissociation and elastic reactions are left as free fit 
parameters. Since only three of the four diffractive cross sections can be reliably extracted from 
the fit we make an additional assumption about the value of the double dissociation cross section 
varying it from to 40 fib. The upper limit is chosen to be about two times larger than the value 



expected from the low energy measurements extrapolated using Regge-type formalism [27|. Our 
attempts to determine the double dissociation contribution directly from the data by observing 
a high mass dissociation of both the proton and the photon give results within this interval, but 
are inconclusive. In order to estimate the contribution of the tail of the non-diffractive reaction 
into the fitted rf m in region, the non-diffractive cross section was fixed in the fit to describe the 
part of the r) m i n spectrum below rj m i n = 1. An example of a fit using the PHOJET model and 
the assumption odd = 20 fib is shown in Fig.||. The fit describes the data well. The cross 
sections a el, &gd arid opd obtained are displayed in Fig.[|a as functions of (Jdd- The errors are 
dominated by systematic uncertainties due to model dependence and are shown by grey bands. 

The model dependence was studied by using different MC generators (PYTHIA and PHO- 
JET) and, in addition, by varying the main parameters of the diffractive model within the MC 
generator. These parameters are the value of the nuclear slope B, the minimum value for the 
mass of the dissociated system and the form of the mass dependence of the cross section: 

• The value of the nuclear slope was varied by AB = ± 4 (GeV/c) -2 in the elastic process and 
half this range in the single dissociation processes. This represents a conservative estimate 
of the uncertainty in the extrapolations of the measured slope from lower energies pi 



The uncertainty in the description of the measured low mass part of the dissociated mass 



spectrum [|l^, 19] was conservatively estimated by increasing the value of the minimum 
mass of the diffractively produced hadronic system by 0.2 GeV/c 2 . 

Mass dependence of single and double diffractive dissociation. The photoproduction data 
at lower energy Jl8| and hadron hadron diffractive dissociation at = 546 GeV ]28| are 
well described by a phenomenological 1/M 2 dependence. This mass dependence is imple- 
mented in the MC models used for the analysis. However, the predictions of Regge theory 
with a supercritical pomeron trajectory give after an integration over t 1/M 2q (°), where 
a(0) is the value of the intercept of the pomeron trajectory at t = 0. The phenomenologi- 
cal fit of the total cross section [Q] gives a(0) = 1.08. In order to estimate this part of the 
model dependence we set the mass distribution also to 1/M 2 - 2 in the MC calculations. 



For every variation of a parameter we average the results and ascribe half of the spread to be 
the corresponding systematic error. This is then added to other errors in quadrature. Similarly, 
we average the results obtained with the PYTHIA and PHOJET models. The final results are 




Figure 5: The measured partial diffractive (a) and non-diffractive and total (b) jp cross 
sections as a function of the assumed double dissociation cross section. The errors are shown as 
bands of ±1 standard deviation. The wide bands correspond to full statistical and systematic 
errors added in quadrature. The narrow grey bands show the systematic errors due to the model 
dependence described in the text. 

shown in FigjBja. The wider bands correspond to full statistical and systematic errors added in 
quadrature. The narrower grey bands show the contribution to the systematic errors from the 
model dependence described above. The main sources of the model independent contribution 
to the systematic errors are the uncertainties in the luminosity measurement (5.6%), in the 
acceptance of the electron tagger (5.0%, which affect all the results) and an uncertainty in the 
statistical background subtraction affecting only the elastic channel (8%). The statistical errors 
are much smaller than the systematic errors in all cases except for proton dissociation where 
they are comparable. Variation of the non-diffractive contribution to the fitted distributions by 
a factor of two, changing the fit interval and variation of the different requirements used for the 
data selection alter the results only within the statistical errors. 

One can see from Fig.|5|a that the elastic cross section is almost independent of any assump- 
tion made about odd- Proton dissociation and especially photon dissociation show a stronger 
correlation with the assumed value of <jdd- However, at any value of ctdd, the single pho- 
ton dissociation is substantially larger than the single proton dissociation, in contrast to the 
assumption that they are equal, made in earlier a^ t analyses (2|, [| at HERA. 

The values of the diffractive cross sections averaged over <tdd> are presented in Table ^. 
8.2 Non-Diffractive and Total Cross Sections 

The measurement of the total 7p and non-diffractive cross sections is based on the ToF-triggered 
sample with the nominal vertex. The efficiencies for the different subprocesses are given in Ta- 
ble |l|. As mentioned above, the efficiencies for elastic and proton diffraction channels using 
nominal vertex are about half those for shifted vertex data making the determination of diffrac- 
tive contributions from the nominal vertex sample less reliable. The diffractive cross sections are 



Table 2: Results for partial and total 7p cross sections under the assumption that the double 
dissociation cross section is in the range < odd < 40 fih. The first error is statistical, the 
second one is systematic, and the third error reflects the systematic uncertainty due to the 
assumption on odd ■ Their sum in quadrature is given as the full error. Note, that due to 
the error correlation in the partial cross sections the error of the total cross section is relatively 
smaller. 



process 




cross 


section (/ib) 






full error (fib) 


o(jp 


-> XY), DD 




20 ± 20 (assumed) 








o(jp 


-> Xp), GD 


23.4 


± 


2.6 


± 


4.3 


± 


10.2 


11.3 


o(jp 


-► VY), PD 


8.7 


± 


1.5 


± 


1.5 


± 


3.0 


3.6 


cr(7p 


-> Vp), EL 


17.1 


± 


1.6 


± 


3.7 


± 


1.4 


4.3 


EL + GD + PD + DD 


69.2 


± 


3.4 


± 


8.8 


± 


9.3 


13.2 




ND 


96.1 


± 


3.5 


± 


14.7 


± 


9.6 


17.9 


Total 


165.3 


± 


2.3 


± 


10.9 


± 


1.3 


11.2 



therefore assumed to be those found from the shifted vertex data, as described in the previous 
section, and the nominal vertex data are used to determine only the missing non-diffractive cross 
section. This was done by solving eq.@ for ond-, then calculating oJ^ t as a sum over five partial 
contributions, all errors being properly propagated including a correlation between diffractive 
and non-diffractive cross sections. 

The <7/vd and oJ^ t obtained as a function of the assumed value of odd are shown in Fig.gb. 
The model dependence displayed has been studied in exactly the same way as in the analysis 
above. The model uncertainty of oJ^ t and ond is dominated by the difference between PHOJET 
and PYTHIA models which enters the calculations via the different efficiencies for the ND, GD 
and DD channels shown in Table |l[ 

The total cross section is remarkably insensitive to assumption about odd and changes only 
by 2.6 fib as odd varies between and 40 fib. This is because the efficiencies e« for GD, DD 
and ND are very similar (see Table [[]). The sum oqd + odd + &nd is therefore practically fixed 
by eq.(||), whilst the total contribution of opd and oel in eq.@ is weakly dependent on odd- 

The results, averaged over odd =0-40 fib are given in Table [2] along with the diffractive 
cross sections. We finally obtain the total photoproduction cross section for an average W 1P of 
200 GeV 

er2S = 165.3 ± 2.3(stat.) ± 10.9(syst.) fib. 

The statistical error reflects all the relevant statistical uncertainties for the shifted and nom- 
inal vertex data samples, as well as those from the MC calculations. The various contributions 
to the systematic errors are listed in Table ||. The dominant sources are the uncertainty of the 
e-tagger acceptance (±8.5/zb) and the difference between the PYTHIA and PHOJET models 
(±5.1iib). 

8.3 Cross Checks of the Total Cross Section Measurement 

The data taken with the Ray-trigger allow a cross check of the model-dependent acceptances 
Ej. Since ToF is sensitive only to particles with 77 < — 2 and the Ray trigger is fired by particles 
around 77 = 0, the two triggers are quite independent of each other. We have compared the 



Table 3: Different contributions to the systematic error of the aj^ t measurement. 



Source 


Error (/A>) 


ToF -trigger 


Ray-trigger 


assumption of odd < 40/ib 

efficiency uncertainty due to model dependence 

syst. errors of diffr. cross sections (without model dependence) 


1.3 
5.6 
1.3 


5.4 
5.7 
4.2 


1.6% uncertainty in the luminosity measurement 

0.4% uncertainty in the fraction of luminosity in satellite bunches 

5% uncertainty in the e-tagger acceptance 

1% uncertainty in QED radiative corrections 


2.8 
0.5 
8.5 
1.7 


Total 


10.9 


12.7 



fraction of ToF-triggered events in the Ray-triggered data sample with predictions from the MC 
simulations. This fraction is 70% in the data, which is 2% higher than the fraction calculated 
with the PHOJET model and 3% lower than that predicted by the PYTHIA model. This cross 
check suggests that the true value of the acceptance lies between the PYTHIA and PHOJET 
estimates, and validates our procedure of averaging the cross sections over the two models. 

A second determination of the cross section has been made using the Ray-trigger data sample, 
giving o~]o t = 162 ± 2 ± 13 tib. Here the relative contributions of the subprocesses are taken from 
the above measurements (section 8.1) for different values of odd- This second determination 
of o~1ot is consistent with the first result described above. The larger systematic error of this 
result reflects the larger difference in the efficiency for different reactions in the Ray-trigger 
data sample, compared to those in the sample taken with the ToF-trigger. The Ray-trigger 
selection has a higher acceptance for the non-diffractive reactions, but shows a nearly vanishing 
efficiency for the elastic and the single proton dissociation processes. In addition, since it depends 
strongly on the transverse momentum and multiplicity of charged particles, the Ray-trigger has 
a higher sensitivity to the details of the hadronic final state simulation. The contributions to 
the systematic error of this measurement are also given in Table ||. 

A determination using only the shifted vertex data resulted in a(^ t = 166 ± 2 ± 15 /ib. The 
larger systematic error is a result of the larger uncertainty of the luminosity measurement. 

9 Discussion 

The energy dependence of a]^ is shown in Fig. ^ with the low energy data |£9[ and the present 
measurement together with a recent result from ZEUS Q . The data are compared with predic- 
tions made by A.Donnachie and P.V.Landshoff (DL) Q and H.Abramovich, E.M. Levin, A.Levy 
and U.Maor (ALLM) ||. The DL curve presents a parameterization of a universal rise of the 
cross section in hadron hadron and low energy photon hadron collisions. In this parameterization 
the high energy cross section behaviour is described by a function with A = 0.0808. The 
dotted line in Fig. || represents a further DL-type parameterization, which takes into account 
the recent CDF measurement of the total pp cross section || with A = 0.11. The ALLM is 
a Regge-type cross section parameterization for real and virtual photon proton collisions with 
A = 0.045. 

In Fig. |H| the measurement of the diffractive cross section contributions is shown as a function 
of the assumed value for the double diffractive dissociation 775 cross section. One should note 




1 10 100 
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Figure 6: Total photoproduction cross section as a function of the 7p centre of mass energy 
Wj P . The solid line is the prediction of the DL [4] combined fit of hadron hadron and low energy 
photoproduction data and the dashed line is the ALLM [5] parameterization. The dotted line 
presents the DL parameterization obtained after the recent measurement of the total pp cross 
section by CDF [8]. 



that the results presented here on photon and proton diffractive dissociation cross sections are 
the first measurements of these quantities at HERA energies. This measurement can be com- 
pared with predictions of A.Capella, A.Kaidalov, C. Merino and J.Tran Thanh Van (CKMT) |30| l, 
G.A.Schuler and T.Sjostrand (SaS) and E.Gotsman, E.M.Levin and U.Maor (GLM) §| 
These models are based on different assumptions about the structure and dissociation of the 
photon and the proton. From these models only the CKMT predictions are based on theoret- 
ical calculations using the Regge model, taking into account absorptive corrections for all the 
diffractive reactions measured here. The comparison of these models with the data is shown 
in Table ^j, where the data are presented for a fixed value of the double dissociation cross sec- 
tion (&dd = 15/xb). This particular choice is made according to the predictions of the models. 
The predictions for the elastic reaction are in good agreement with the HI measurement as well 
as with the recent result from the ZEUS collaboration (18 ± 7/xb) |J. However, the relative 
contribution of single proton dissociation is observed to be about three times lower than that 
of photon dissociation. This observation disagrees with the predictions of the SaS and GLM 
models. In the latter the ratio of the proton to the photon single diffractive dissociation is 
obtained using the quark counting rule. The results of the CKMT calculations are supported 
by the present measurement. 



The published results of the GLM calculation are scaled to a dissociation mass interval of M 2 < 0.1 



Table 4: The comparison of the diffractive cross section calculations with the HI measurement. 
The data are presented for a fixed value ctdd = 15/zb. 



Reaction 


Cross Sections (/ib) 


Data 


CKMT 


SaS 


GLM 


a (HP ~^ Vp), EL 


17 ± 4 


17 


16 


17 


<r(lP^ Xp), GD 


26 ± 5 


25 


13 


18 


o-(7P-» VY), PD 


9± 2 


7 


10 


15 


a(-f P ^ XY), DD 


15 


15 


13 


15 



10 Conclusion 

Using the HI detector at HERA results on 7p scattering at the average cm. energy of W 7P = 
200 GeV are obtained. The total photoproduction cross section is measured to be oy/ot = 
165 db 2 (stat.) ± 11 (syst.) fib replacing our previous result |^]. The extracted diffractive 
and non-diffractive contributions to the cross section are presented as a function of the assumed 
value for the double diffractive dissociation cross section odd • The cross section of single photon 
diffractive dissociation is observed to be substantially higher than the single proton diffractive 
dissociation cross section. The elastic cross section depends only weakly on the assumption 
made on odd and is found to be 17.1 ± 4.3 ^b . 
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